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Due to their collective nature Casimir forces can strongly depend on the geometrical shape of 
the interacting objects. We study the effect of strong periodic shape deformations of two ideal 
metal plates on their quantum interaction. A non-perturbative approach which is based on a path 
integral quantization of the electromagnetic field is presented in detail. Using this approach, we 
compute the force for the specific case of a flat plate and a plate with a rectangular corrugation. 
We obtain complementary analytical and numerical results which allow us to identify two different 
scaling regimes for the force as a function of the mean plate distance, corrugation amplitude and 
wave length. Qualitative distinctions between transversal electric and magnetic modes are revealed. 
Our results demonstrate the importance of a careful consideration of the non-additivity of Casimir 
forces, especially in strongly non-planar geometries. Non-perturbative effects due to surface edges 
are found. Strong deviations from the commonly used proximity force approximation emerge over a 
wide range of corrugation wave lengths, even though the surface is composed only of flat segments. 
We compare our results to that of a perturbative approach and a classical optics approximation. 

PACS numbers: 03.70.+k, ll.10.-z, 42.50.Ct, 12.20.-m 



I. INTRODUCTION 

Casimir interactions Q, H II H are a fundamental property of the vacuum. They are commonly related to 
quantum electrodynamics but fluctuation induced interactions are of interest in a plethora of other fields like in 
condensed matter systems as liquid crystals and superfluids 0, 0> m cosmological models ||, in particle physics 
d llf)| | and in biological systems as proteins on membranes. The guiding mechanism behind all these phenomena is 
that a quantum or thermal field is constraint by boundary conditions on surfaces so that the energy is modified, and 
effective interactions between the surfaces occur. For quantum fields the Casimir interaction is given by the change in 
the ground state energy £o — i h\o n due to the presence of boundary conditions. Even for non-interacting fields 
like the photon gauge field it is difficult to obtain the Casimir interaction since the eigenfrequencies uj n can depend 
strongly on the confining geometry. Thus it is not unexpected that exact analytical results for Casimir interactions 
between macroscopic objects are not known even if the geometry has high symmetry. 

Most of the recent high precision experiments aim at the measurement of the Casimir force in geometries which are 
closely related to the standard case of two parallel plates [H El El El To the latter case applies Casimir's seminal 
prediction [l| 
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for the force between two ideal metallic and parallel plane plates of area A and distance H at zero temperature. For 
technical reasons, usually a plate-sphere geometry is used in experiments. Even for this case only approximative 
methods like the proximity force or Derjaguin approximation [15j can be applied and the exact result is not known. 
There exists geometries for which there is even little intuition as to whether the interaction is attractive or repulsive. 
A striking example is Boyer's result that the Casimir energy of a conducting sphere is positive |Tfij . This observation 
has triggered a search for repulsive configurations [l?! Hsj . Such effects can be even of direct practical relevance in 
nano technology where "sticking" of mobile components in micromachines might be caused by Casimir forces |f 9j. 

The advances in experimental techniques have stimulated the measurement of the shape dependence of Casimir 
forces in specially designed geometries (as opposed to inevitable geometrical effects such as surface roughness). Mo- 
hideen et al. were able to measure the Casimir force between a sphere of large curvature radius and a corrugated plate 
[20L l2lj . Although the corrugation length was larger than the studied range of separations between the surfaces, their 
results showed a clear deviation from predictions of the proximity force approximation. While it has been suggested 
[22l | that lateral surface displacements caused this discrepancy, there is no reason to believe in the validity of the 
proximity approximation if the corrugation length is decreased. 

Because of the wide range of realizations of Casimir forces, improved experimental techniques and the increasing 
importance of nanostructures it is interesting to develop novel approaches for computing such interactions. In the limit 
of slight surface deformations a p ath integral quantization subject to boundary conditions allows for a perturbative 
calculation of the interaction pa |24| , showing strong corrections to the proximity approximation |29 | . Another 
perturbative approach, based on a multiple scattering expansion, has been applied to the limit of large surface 
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separations |2jfl . Very recently, an alternative approximation scheme based on geometric optics has been proposed for 
geometries where the Casimir interaction is mostly caused by short wavelengths |2^. However, to date no systematic 
method is known for estimating interactions of strongly deformed objects, including large curvature or even sharp 
edges. In this paper we present a novel non-perturbative method to compute electrodynamic Casimir interactions 
between uniaxially and periodically deformed surfaces. It is based on a path integral approach for Casimir forces 
[2fj, l3fT | . The approach is not restricted to small deformations or small surface curvature but it also allows us to 
study strong deformations and edges. We develop a numerical implementation of the approach which allows for a 
precise computation of the interaction without any approximations. As an example we consider a geometry consisting 
of a flat and a rectangular corrugated plate, see Fig[21 For this geometry we obtain the Casimir force over a wide 
range of surface separations and corrugation lengths. We find that the edges of the corrugated surface cause strong 
deviations from the proximity approximation which agrees reasonably with our results only if the corrugation length 
is much larger than both the surface distance and the corrugation amplitude. We show that the qualitative effect of 
edges on the interaction can be understood in the limit of large corrugation lengths in terms of classical ray optics. 
A brief account of our method and its application to scalar fields subject to Dirichlet boundary conditions appeared 
in Ref. EH 

The rest of this work is organized as follows. In Sec. II we review briefly the path integral approach and then 
introduce the method for a non-perturbative computation of Casimir interactions. We consider periodic uniaxially 
deformed surfaces. The new approach is then applied in Sec. Ill to the example of a flat and a corrugated surface 
with sharp edges. The asymptotic limits of small and large corrugation lengths are treated analytically. For arbitrary 
corrugation lengths the interaction is obtain by an numerical implementation of our approach. We give detailed 
numerical results for the total electromagnetic Casimir force and the contributions from TM and TE modes separately. 
In Sec. IV we compare our results to perturbation theory for slightly deformed smooth surfaces. We interpret our 
results for large corrugation length in terms of classical ray optics. Throughout the paper we set c = 1 and H = 1. 



II. NON-PERTURBATIVE PATH INTEGRAL APPROACH 



We consider two perfectly conducting periodically deformed (corrugated) plates S a , (a = 1,2) with a mean sepa- 
ration H. They are assumed to be infinitely extended over a base plane which is parameterized by the coordinates 
X|| = (xi, £2). For simplicity we assume that the corrugation is uniaxial along the x\ direction. The shape of the plates 
is then described by height functions h a (x\) which measure deviations from the mean height so that J h a (x\) = 0. 
The Casimir energy of the two plate configuration can be obtained from an imaginary time path integral representation 
[23T,l24| for the partition function of the electromagnetic field and the confining plates. In the absence of boundaries, 
the path integral extends over the electromagnetic gauge field with the 4D space time action 

So [A"] = 7 / d 4 XF^F^. (2) 
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and the field — d^A v — d u A^ and X = (xq, X|| , X3). In order to eliminate redundant gauge field configurations 
the Fadeev Popov gauge fixing procedure has to be applied [12 ■ The ideal metal boundary condition for the gauge 
field A^{X) is given by the requirement that the tangential components of the electric field vanishes at the surfaces. 

For plate deformations which are uniaxial, the translational invariant direction can serve as reference axis for defi ning 
transverse magnetic (TM) and transverse electric (TE) modes, similar to the treatment of wave guide geometries [33j . 
Then every field configuration can be decomposed into those two types of modes, and one can resort to a scalar field 
path integral quantization 25, 26|. The scalar fields are given by the electric and magnetic field component along the 
translational symmetry axis, 

$(A) = E 2 {X) for TM modes, (3) 
$(A) = B 2 {X) for TE modes. (4) 

Since the plates are assumed to be ideally conducting, the boundary conditions for TM and TE modes are of 
Dirichlet and Neumann type, respectively, i.e., 

$| s =0 for TM modes, (5) 

<9 fi J>| s = for TE modes, (6) 

with the surface normal derivative denoted by <9„ pointing into the vacuum between the plates. After a Wick rotation 
to imaginary time x$ — > ix$, both types of modes are described by the Euclidean action 

Se{*} = \ J d 4 A(V$) 2 . (7) 
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In 4D Euclidean space, the surface positions of the plates are then parameterized by X a (r) — [r, h a (xi) + H8 a 2\ with 
r = (xo,Xii). Following the procedure introduced in Refs. |2iJ |3(], the boundary conditions are imposed by inserting 
delta functions on the surface in the functional integral. The partition function for TM and TE modes, respectively, 
then reads 

Zd = [ ®<S> lH[5l$(X a )] e- s *«, (8) 

J a=l X a 

Zn = z^ 1 f m J] Y[5[d & j>(x*)} ^ SeW (9) 



a=l X a 



with the boundary free partition function Zq. The functional integrals can be calculated by introducing auxiliary 
fields to represent the delta functions. Then, the Gaussian integration over $ can be carried out, yielding the partition 
function in terms of an effective action for the auxiliary fields, 



r 2 



_fie- 4 " W °' (10) 

a=l 

with the effective action 



Seff{><*} = \ J /£ Mr)M a "(r; r')#(r'). (11) 

In the following we will drop the subscript D or N for the boundary conditions on all quantities which apply to both 
conditions in the same way. The total electrodynamic Casimir energy is then given by the sum of TM and TE mode 
contributions, E = £tm +£te- After subtracting the divergent and H independent terms, the energies can be written 
as £tm = lndet(A / tDA / tp 1 oc )/(2AL), and analogous for £te with D replaced by N, where Air>,oo is the asymptotic 
expression of A^d for H — ► oo, A is the surface area of the plates, and L is the Euclidean length in time direction. 
The Casimir force F = —dn£ per unit area is then given by F = -Ftm + Fte with 

Ftm = -^-^(M^dHMu), (12a) 

Fte = -^^W^n). (12b) 

The right hand side of these expressions is always finite, and no regulation of divergences by subtraction of the 
vacuum energy in the absence of boundaries is necessary. The Dirichlet and Neumann matrix kernels of the effective 
Gaussian action can be expressed in terms of the Euclidean scalar Green function Q(r, X3) = (r 2 + x 2 )^ 1 j (47r 2 ), and 
are respectively given by 

M^(r;r') = Va ( Xl ) m (x[) Q (X a (r) - X (r')) , (13a) 
(r; r') = Va ( Xl ) Vp {x{) d^^d^^Q (X a (v) - X p (v')) (13b) 

with the coefficients given by r) a (xi) = (l + [(d Xl h a )(xi)] 2 ) 4 . These coefficient arise from the integral measure on 
the curved surfaces. However, since they are independent of the mean plate distance H, they cancel in the matrix 
product of Eq. (|12fl and therefore can be ignored for the calculation of forces. The kernels are symmetric with 
A4(r; r') = Ai T (r'; r) where the transpose refers to a, (3. Using the parameterization in terms of height profiles, the 
matrix kernels can now be written as 

(r; r') = Q (r - r', h a {x x ) - hp(x{) + H(5 a2 - 8 fj2 )) , (14a) 
A<(r;r') = (-1)^ [-d 2 X3 + (h' a ( Xl ) + h'^x',)) d Xl d X3 - h'Mhfai) d 2 Xl ] 

x g (r - r', x 3 - 4) X3 =h a (xi)+HS a2 ( 14b ) 

x' 3 =h fj (x' 1 )+HSp 2 

for the Dirichlet and Neumann case, respectively. So far, we have not used the periodicity of the surface profile, and 
the above results are valid for any uniaxial deformation. However, the computation of the force can be performed 
more efficiently if the periodic symmetry of the surface is used. 
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FIG. 1: (color online) Transformation of the matrix M to block-diagonal form. The figure shows a finite part of the matrix, 
corresponding to the blocks Bki with k, I — —1,0,1, before and after the permutations of rows and columns. Before the 
transformation (left box) M has a band structure with diagonal blocks Bij consisting of 2 x 2 matrices N m along the diagonal. 
(The dependence on the lateral momentum qx is not shown here.) The first step of the transformation is to permute the 
rows and columns which are formed by the first entry N m in every block Bki (indicated as grid). These entries from after the 
permutations the first block Mo of M (right box). The latter permutation process is then repeated for the second entry, the 
third entry until the (N + l) th entry of every block Bki, leading to the N + 1 blocks Mj. The momenta qi within each block 
Mj are constant for every column and they differ only by integer multiples of 2n/X between columns (of the same block), see 
labels in the right box. The blocks Mj differ in their momentum shift jS, 5 — 2ir/W , which is located in the unit cell [0, 2n/X[ 
since j = 0, ... ,N = W/X - 1. 



Due to the translational invariance in time (xq) and one space (2:2) direction, it is convenient to introduce the 
momentum vector qj_ = (90 5^2) which is perpendicular to the direction of modulation. Due to the periodicity of the 
surface profile, the Fourier transform M (p; q) = f f , e i P r + M i r M.(r\ r') can be decomposed into the series 

00 

M(pi,px;gi,q±) = (27r) 3 (5(p ± +q±) ^2 S (p x + q x + 2nm/X) N m (q x , q x ) (15) 

m— — 00 

where N m (q±,qi) arc 2x2 matrices which depend only on q± = |q_i_|. From Eq. (|15fl it is obvious that the matrix 
M. has its non-zero entries arranged in 2 x 2 blocks along parallel bands. Due to this structure, there exists a 
transformation, consisting only of row and column permutations, which makes the matrix block diagonal. To perform 
this transformation, we cut the matrix M. into smaller matrices B k i which have non-zero entries only in 2 x 2 blocks 
along the diagonal, see Fig^ For the purpose of parameterization, we consider discrete momenta p x = (2ir/W)j, 
j = 0, . . . , N, along the direction of surface modulation with N — W/X — 1. The continuum limit is obtained if the 
linear size W of the surfaces and N are taken to infinity in order to obtain the force per surface area A = W 2 . With 
this parameterization, the block diagonal matrices B k i of dimension 2(N + 1) x 2(N + 1) can be read off from Eq. JTSJ, 
leading to 

B k i(q ± ) = dia,g{B kl (q ± ,0),B kl (q ± ,2Tr/W),...,B kl (q ± ,2nN/W)} (16) 

with the 2x2 block matrices defined as [see Fig. Q] 

B kl (q ± , qi ) = N k _i (q ± , qi + 2-kI/X) . (17) 

By inspection of Fig.n° nc easily realizes that a sequence of row permutations and a subsequent sequence of column 
permutations transforms the matrix M. to block-diagonal form. Each of the N + 1 blocks M.j is composed of exactly 
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one element from each matrix Bki and those elements forming a block Mj come from the same position in every 
matrix Bki as indicated by the color scheme in Fig.^ Thus, each block Mj is composed of entries which correspond 
to the same discrete momentum pi = (2ir/W)j, and we obtain for the elements of Mj the result 

MjM (P-L><U-) = (2tt) 2 5 (p ± + qx) B M (q ± ,2irj/W) . (18) 
The number of permutations needed for the matrix transformation is always even, and thus we get the determinant 

JV 

detM = JJdetAlf. (19) 

3=0 

By differentiating with respect to the mean surface distance H and by using the relation lndetMj = TrlnMj, we 
obtain 

N 

d H (in detM) = J^Tr (Mj 1 ■ d H M^j . (20) 

3=0 

This result reflects the fact that the (free) energy of the system can be calculated as the sum of the individual (free) 
energies of decoupled subsystems, which are described by the matrices Mj. Each subsystem with fixed j describes 
scattering events at the fixed momenta pi = (2ir/W)j + (2tt/X)1 which differ only by integer multiples of 2tt/\. 

Using Eq. (|18f) we can perform the trace over the continuous lateral momenta and the discrete indices within a 
fixed subsystem, 



Tr (Mj 1 ■ d H Mj) = J d 2 ^ 1 - B M, a p (?-L, Qi) ' d H B lk j ja (q±,qi) 



(21) 



k,l— — oo 
a, (3=1,2 



where we have explicitly indicated that the trace is performed with respect to all discrete indices, and we remind that 
L is the system size in time direction. It appears useful to define the function 

g (q ±>qi ) ee tr (B'^qx, qi) ■ H B (q±, qi )) , (22) 

with the lower-case symbol tr denoting the trace over the discrete indices summed over in Eq. I|21() . Next we perform 
the sum over all subsystems with j = 0, . . . , N — W/X — 1. This can be easily done by going back to continuous 
momenta p±. If we take the limit W, N — > oo with W/ (N + 1) = A fixed the sum in Eq. I|2()(l can be written as the 
integral 

d H (In detM) = -— ^ / d 2 q±— / dq ig (q ± , qi ). (23) 



(2^) 2 J ^ 2n J 

The function g(q±, qi) has the following symmetry properties. A shift of the momentum p\ by 2tt /X corresponds just 
to a renumbering of the matrix elements Bki since the matrix is of infinite dimension. Thus we have g(q± , qx + 2tt/X) = 
g(q±, <7i). If both surface profiles are described by even functions, h a (— Xi) — h a (x\), for the matrices N m the relation 
N m (q±, —qi) = N_ m (q±,qi) holds. Using the later relation and the definition of B^j of Eq. I|17fl it is easy to check 
that g(q±,—qi) = g(q±,qi) by performing appropriate row and column permutations for the matrix B. The above 
symmetries allow to write the Casimir force per unit area, F/A — — (\/2LW 2 ) 3//(lndetM), as 



1 /"OO />7r/A 

F l A = -^J n d( l^ j dq ig ( q± , qi ). (24) 



This is the final result of the general approach for arbitrary uniaxially corrugated surfaces. As we will show below, it 
can be used for an efficient numerical computation of the Casimir force. The input of such a numerical approach are 
the matrices N m from the decomposition in Eq. I|15l) . Moreover, the result can be also used to obtain non-perturbative 
analytical results in the asymptotic limit of very small corrugation lengths. 

Before one can develop a numerical implementation of the above representation of the Casimir force, one of course 
has to restrict the infinite dimensional matrices. In the remaining part of this section we will introduce a suitable cutoff 
procedure for the matrix dimension. We will take two flat plates as a simple example to examine the convergence 
of the procedure if the cutoff is taken to large values. The cutoff procedure consists in the restriction of the matrix 
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M to blocks B M with ft, I = —M, ...,M only. The dimension of M is then 2(2M + 1)(N + 1). Fig. [T] displays the 
restricted matrix for M = 1. The corresponding function g is then defined by Eq. 12 I'll with the restriction that the 
trace runs over ft, I = —M, . . . ,M only. We will denote this function in the following by gu- This function is then 
used instead of g in Eq. (|24[l to obtain a series of approximations Fm to the force which converges to F for M — ► oo. 
As an example consider two flat plates at distance H . Then M. is a diagonal matrix and N m = for m 7^ 0. Thus 
the matrix i? is also diagonal with 

Bkl(qx,qi)=SklN (q ± , qi +2nl/\). (25) 
Using Eq. (1221) with the trace taken for ft, Z = —M, . . . , M, one gets the function 

^ 2y/ g i + fa + 27Tl/\r 

9Miq± > qi) = j ^i-H^A)^-! ' (26) 
Integration over gi yields an M th order approximation Fm to the force, 



1 



2tt/A 



I r°° /-27r(M+l)/A 2 g 

8?H Jo J-2kM/\ z q - 1 



with q = v/3j_ + For — > 00 one gets the known (A independent) result F/A = —(tt 2 /480)H 4 , and the finite 
M corrections to this asymptotic result scale exponentially fast (~ e -^MH/\^ ^ Q zero f or i ar g e y[_ Therefore, in the 
case of periodically deformed plates one can expect accurate numerical results for F from moderate values for the 
cutoff M, and the convergence is faster for smaller A. 



III. RECTANGULAR CORRUGATION 



In the previous section we developed a non-perturbative approach for computing Casimir interactions between 
periodically deformed surfaces. In this section we will use the approach to obtain explicit results for the Casimir force 
between a flat plate and a plate with a rectangular grating. The effect of this class of periodic geometries (corrugated 
surfaces) can significantly modify the interaction of the objects |34l l35j . Is was proposed that such geometries can be 
used to reveal novel features of the Casimir interaction |23l l24j . For a similar geometry consisting of a sinusoidally 
corrugated plate and a sphere with a radius S> H Roy and Mohideen measured the force, and found clear deviations 
from the predictions of the proximity force approximation 20] . While it has been suggested that lateral shifts of the 
surfaces caused the discrepancy, we demonstrate below that periodic surfaces allow for a much stronger sensitivity 
to geometry if the corrugation length is reduced to smaller values. Specifically, we consider the geometry shown in 
Fig. [5] with a rectangular grating of amplitude a and wavelength A. Choosing x\ as the direction of modulation, this 
corresponds to the height profile 

= { +a J " 'wJ < ^1 wo. (28) 
iV 11 \ -a for A/4 < |sei| < A/2 ' v ; 

and continuation by periodicity hi(xi) — hi(xi + nX) for any integer n. The upper plate is flat so that /^(^l) = 0. 

The main purpose of our work is to obtain the Casimir interaction in regimes where other methods like proximity 
approximation, pairwise summation of two body forces or perturbation theory fail or become unreliable. While the 
proximity approximation assumes smooth profiles with small local curvature also perturbation theory in the height 
profile yields divergences in the presence of edges in the profile or if the corrugation length becomes very small, i.e., 
A <C a, H [25|,|2(|. In perturbation theory, finite corrections of order a 2 are recovered only if edges are "smeared out" 
over a finite length scale. Thus the correct procedure would be presumably to sum all orders of perturbation theory 
for a "smeared out" profile, and then to take the limit of sharp edges after summing all contributions. Since the 
perturbative treatment is rather cumbersome the latter program is not practicable, and non-perturbative techniques 
are imperative. 

In order to apply the general result of Eq. H24|) for the Casimir force, we have to decompose the matrix M. into the 
matrices N m , see Eq. I|15fl. For a general profile, this has to be done by a numerical Fourier transformation. A nice 
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FIG. 2: Geometry consisting of a rectangular corrugated plate and a flat plate. The surfaces are translationally invariant 
along the X2 direction. 



property of the rectangular profile of Eq. I|28[) is that it allows for an analytical computation of the matrices N m . The 
idea is to rewrite the profile of the corrugated plate as a discrete Fourier series, 



2 °° 
hi(xi) = — V 

7T ^ ' 



1 c 2fi(2n-l)xi 



7T <—< 2n - 1 

n— — oo 



(29) 



which is inserted into the matrix Ai of Eq. (|14() . Then the Fourier transformed matrix Ai can be calculated, leading 
after some algebra to the matrices N m (q± , q\ ) . Details of this calculation are given in Appendix A. The results are 



A%(q±,qi) 
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I V 

for Dirichlet conditions and 

' / Af n (q ± , qi ) 

V 

Nji tTn (q±,qi) = < 



(-1)" 







2q 

(-1)^ e-" H 



+ 0-mO 



m.7T q 
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sh(ag) 
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for rn even 



for m odd 



(30) 



f(l + e- 2a 9) §e^ H ch(a 9 ) \ 



%e-i H ch(ag) 



(-1) 2 r -q m .H 2ttto qi+27rm/A -| 



m A 



for Neumann conditions with 



^m(<U,<7l) = "2 



sh(ag m ) 



-2aij 2 fc-l _ ^ 



7T- , ^ (2fc- l)(m-2fc + l) g 2fe _i 

/c— — oo 



for m even 

for m odd 
(31) 

(32) 



and 



An (<7-L,<7l) 



E 



M) 1 



fe=-c 



X _ g-2aij2fc 



(2fe-l)(m-2fc + l) g 2 Vr 



(33) 



/ 27rm , , 27r . , . , 2 „ 9 / 7rm N , 27T . s » 

Qi (qi + — ) (91 + — (2* - 1)) + 2q 2 ± ( Ql + — ) ( 9l + — (2k - 1)) + q! 
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respectively, with the definition q n = \Jq\ + {qi + 2Tm/X) 2 which implies q = qo- With these results at hand, the 
Casimir force can be calculated by the approach developed in the previous section. The recipe is as follows. First, 
one constructs the matrix Bki of Eq. I|17|l. then one calculates the inverse of Bki to obtain the function g(p±,pi) of 
Eq. (|22[) . and finally one has to perform the integration of Eq. (|24[) . In general, this program can only be performed 
numerically. However, in the limit A — > it turns out that a closed form for the function g(pj_,p) is available which 
allows to obtain the Casimir force in this limit exactly. 



A. The limit of small A 



Let us consider the case where the corrugation length A sets the smallest length scale in the geometry of Fig 13 If 
we take the extreme limit of A — > 0, a naive assumption is that the field can no longer get into the narrow valleys 
of the corrugated plate. Even for small but finite A this picture should be a good, though approximate, description 
since it still effects the wavelengths of order H which give the main contribution to the force. Thus one expects that 
the plates feel a force which is equal to the force between two flat plates at the reduced distance H — a. However, the 
question remains to what extent this is a good approximation when A becomes larger, say of order a. To check our 
naive expectation, we will apply the approach of the previous section to the limit A — > 0. Fortunately, in this limit the 
matrices N m (q±,qi) simplify considerably both for TM and TE modes. The explicit form of these matrices is given 
in appendix B. From this result, we can explicitly calculate the functions <7m(9_l, <7i) which was introduced before 
Eq. (|25|) . As explained in section [HI the infinite dimensional matrix Bki is truncated for the calculation at order 
M with k, I = —M, . . . , M so that the truncation is done symmetrically around the center at (k, I) = (0, 0) which 
contains the leading matrix entries. From the exponential convergence behavior of the flat plate result given below 
Eq. (|2Tjl one can expect that in the extreme limit A — *■ 0, the series gM(q±,qi) converges so rapidly towards g(q±,qi) 
that already for M = 1 the exact asymptotic expression is obtained. Indeed, our explicit calculation of <?m(q_l,9i) for 
low M confirms this expectation. From the truncated matrix Bki of Eq. I|17|) and the matrices of appendix B we get 
the simple result 



-2a<n 



i ' 1 , r 2.,_, e 2(H-. ) , for M — 
9M{q±,qi) = < r 1+ \ , (34) 

{ q [coth {q{H - a)) - 1] for M > 1 

for both TM and TE modes. Thus from first order (M = 1) on the function <7a/(<7j_; <7i) remains invariant with 
increasing dimension M of the matrix Bki- Interestingly, the result for M > 1 has precisely the form, which one gets 
for two flat plates at reduced distance H — a. If one integrates the function <7m(<7_l, <7i) for M = 1 one obtains from 
Eq. I|24|) the Casimir force per surface area 

F ° /A = -WoW^W (35) 

for both TM and TE modes. Thus in the limit A — > both types of modes yield the same contribution to the total 
electrodynamic Casimir force F = 2Fq. The result of Eq. (|35fl corresponds to the naive reduced distance argument 
given at the beginning of this section. Notice that this result is non-perturbative in a/H and is exact in the limit 
A — > 0. Perturbation theory for smoothly deformed surfaces always yields corrections to the force of order a 2 [25ll26| . 
However, for small a/H, the result of Eq. (|35|l has the expansion 



F °I A - "on m 



480 fl 1 



4— + 0[ ( — 

H \\H 



(36) 



which indicates that perturbation theory is not applicable if A -C a. Below we will see that the force Fq provides an 
upper bound for the Casimir force from both TM and TE modes at fixed H/a, i.e., for increasing A the force always 
decreases compared to Fo . We expect that the results of this section for A — > are valid for corrugations of arbitrary 
shape and also for rough surfaces if A is identified with the characteristic length scale for surface deformations. 



B. The limit of large A 



In the opposite limit of very large A the corrugated surface is composed of large flat segments with a low density 
of edges. At sufficiently small surface separations H <C A the main contribution to the force comes from wavelengths 
which are much smaller than the scale A of the surface structure. Thus in the dominant range of modes diffraction 
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can be neglected, and the simple proximity force approximation (Derjaguin approximation [l^]) should be applicable. 
Such an approximation assumes that the total force can be calculated as the sum of local forces between opposite 
flat and parallel small surface elements at their local distance H — h(x\). No distinction is made between TM and 
TE modes. This procedure is rather simple for the rectangular corrugation considered here since the surface has no 
curvature (except for edges). There are only two different distances H + a, H — a which contribute one half each 
across the entire surface area, leading for A — * 00 to the proximity approximation for the force, 



Foo/A 



ir z 1 
480 2 



1 



1 



(H - a) 4 (H + a) 4 



(37) 



Below we will see that later result provides a lower bound for the Casimir force from both TM and TE modes, 
contrast to the limit of small A the correction for small a/H is of order (a/H) 2 here. 



In 



C. Numerical analysis 




FIG. 3: (color online) Casimir force for TM modes (a) and TE modes (b) as function of H/a for different corrugation lengths 
A/a. Displayed is the change of the force compared to the force between two flat plates, ^TM.flat = -FTE.flat = 
in units of i*TM,flat and -FxE.flat, respectively. The two bold curves enclosing the numerical data are the analytical results Fo 
for A — > (upper curve) and Fx, for A — > 00 (lower curve), see text. 



In this section we implement the non-perturbative approach of sectionOnumericallv for the rectangular corrugation 
of Fig. [3 One has to resort to a numerical analysis here since the function g(q±, qi) cannot be obtained analytically 
from the matrices of Eqs. (|30[1 . i|31|) for arbitrary corrugation lengths A. The numerical procedure follows straightfor- 
wardly the computation of the Casimir force in section [H] The following implementation applies both to TM and TE 
modes. At fixed order M, the truncated matrix Bki of Eq. I|I7|) with k, I — —M, . . . , M is calculated from the matrices 
N m of Eqs. (|30[1 . Then the matrix B^i is inverted numerically to yield the function <?m(<Z_l,9i) from Eq. H22|) 

where the index M denotes the truncation order. Notice that the derivative of Bki with respect to H is obtained 
analytically and no potentially inaccurate numerical derivatives have to be computed. Finally, the integration in 
equation (|24() is carried out numerically without difficulty since <?m (<Z_l 1 <Zi ) decays exponentially fast for large q± , qi . 
This provides a series of approximations Fm to the Casimir force which must converge to the exact value of the force 
as M — > 00. From our analysis of the flat plate geometry, see Ea. (|27|l . we expect an exponentially fast convergence 
F — Fm ~ e~ lM with a coefficient 7. However, the decay coefficient 7 depends on the geometrical lengths, and it 
is expected to increase with decreasing X/H. This type of convergence behavior we found to be consistent with our 
numerical data for Fm- It allowed us to extrapolate the data to obtain the Casimir force F. The largest M for which 
we calculated Fm was M = 10 for small A/a = 0.1 and M = 97 for large A/a = 300. 
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X/a = OA 
X/a = 1 
X/a =5 
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XI a = 30 
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FIG. 4: (color online) Total Casimir force as sum of TM and TE mode contributions in the short distance regime. Shown is the 
relative change of the force compared to the total Casimir force Ffl at between two flat plates. The data enclosing bold curves 
have the same meaning as in Fig. [3] but for the total force they are now given by 2Fo and 2Foo due to the same contribution 
of TM and TE modes in these two limits. 




FIG. 5: (color online) (a) Ratio of Casimir force from TM and TE modes as function of the plate distance H for different 
corrugation lengths A. (b) Logarithmic plot of the deviation of the ratio from one at large H. 



The results of our numerical analysis are as follows. If we express the total Casimir force F or the force contributions 
-Ftm and Fte from TM and TE modes, respectively, in units of the corresponding force between two flat plates the 
results can be expressed in terms of the dimensionless ratios H/a and A/a only. The results from the extrapolation 
of the data for Fm are shown in FigOboth for TM and TE modes and different corrugation lengths. For both types 
of modes the force Ftm , -Fte is bounded at a fixed plate separation H/a between Foo and Fq as given by Eq. IpTfjl 
and Eq. respectively. For small A/a the upper bound Fq is approached whereas for asymptotically large X/a 

the force converges towards the lower bound F^ which is given by the proximity force approximation. Since the 
convergence towards the lower bound i 7 ^ becomes slower with increasing H/ a there are two distinct scaling regimes 
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FIG. 6: (color online) Scaling of the force from TM (a) and TE (b) modes close to the upper bound Fo (A — » 0) and the lower 
bound Foo (A — > oo) as a function of A/a at fixed mean surface distance H = 10a. Forces are measured in units of FrM.flat and 
-FTE.flat, respectively. 
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FIG. 7: (color online) Same plot as in Fig. but for fixed distance H = 100a. 



for the force at a fixed corrugation length X/a. At small H/a the relative change of the force compared to the force 
between two flat plates, i r T/-fT,flat — 1, T=TM or TE, decays as (H/a)~ 2 . After a crossover regime the relative change 
of the force decays at larger H ^> A like (if/a) -1 , following the behavior of the exact result Fo for A — > 0. The so far 
described qualitative behavior of the force is common to both types of modes. However, there is a clear distinction 
between TM and TE modes, especially at large A/a, as can be seen from Fig|3J The force from TE modes has much 
more pronounced deviation from the proximity approximation result Foo as the TM modes. In particular at large 
corrugation lengths (A/a = 300) this can be seen clearly from our numerical data. The same behavior is observed 
for the deviations from Fo at small A/a. Thus, the force Fte appears at intermediate values of A/a more strongly 
separated from the lower and upper bound, cf. Fig- El^b). We will come back to this point below when we discuss the 
scaling of the force with A close to the bounds. Fig. ^ shows the total Casimir force in the range of small separations 
H. 
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For particular geometries like a cubic volume the Casimir force has even a different sign for a scalar field with 
Dirichlet boundary conditions (attractive force) and an electrodynamic field (repulsive force) 0, |3(|. Since for 
uniaxial plate deformations both types of fields differ in the presence of a scalar field with Neumann boundary 
conditions (TE modes) it is interesting to study more quantitatively the difference between the two wave types. Fig- El 
shows the ratio Ftm /Fte of the forces from both types of modes at different A/a. One observes that the ratio is peaked 
at a characteristic H/a which depends on A/a. For small H/a — ► 1 the ratio tends to one as one can expect from 
the proximity force approximation which does not differentiate between TM and TE modes. In the opposite range 
of large H/a again both types of modes must contribute almost equally since the geometry approaches that of two 
flat plates. For the entire range of studied corrugation lengths the ratio converges to one for large H/a according to 
\Ftu/Fte — 1| ~ (H/a)^ 1 , see Fig.|5Jb). However, this asymptotic behavior sets in only beyond a crossover separation 
H which increases with A. At intermediate A/a the ratio varies approximately between 0.95 and 1.15 in the studied 
range of A/a. TM modes dominate at A/a < 10 and at small H/a for all A/a. The contribution from TE waves is 
larger for A/a > 10 and H/a > 2. It is instructive to compare this behavior to perturbative results of Ref. EME^I for 
the geometry consisting of a smooth sinusoidally corrugated and a flat plate. As will be explained in more detail in 
the next section, the perturbative result for the later geometry yields Ftm/Fte > 1 for all A/a 3> 1 and H/a ^> 1, in 
contrast to our results for the rectangular corrugation. This observation suggests that the corners of the rectangular 
corrugation in fact cause the slight amplification of TE modes compared to TM waves at A/a > 10. One can argue 
that imposing for TE modes a vanishing normal derivative on the field at the concave corners inside the valleys of the 
corrugation provides a stronger constraint on field fluctuations as compared to Dirichlet conditions for TM modes. If 
the width of the valleys is decreased with A the two opposite corners can no longer be considered separately and the 
Dirichlet condition might provide a stronger restriction. For very small H/a the main contribution to the force comes 
from rather short wavelengths which should be only very weakly effected by the Neumann conditions at the concave 
corners. 

Finally, we consider the scaling of the force from TM and TE modes close to lower and upper bounds F^, and 
Fo, respectively. Figures El and show a logarithmic plot of force form TM and TE modes at fixed H = 10a and 
H = 100a, measured relative to for large A/a and relative to Fo for small A/a. At small A we found an interesting 
qualitative difference between TM and TE modes for the scaling towards the exact result Fo for A — * 0, 

F — Ftm ^ A Fq — Fte ^ / A \ 12 
Ftm, flat a' Fte, flat \ a / 

For the change in the exponents we cannot present a satisfying simple argument. In the opposite limit of large A the 
proximity approximation result Foo is approached linearly for both types of modes, 

Ftm — Fqo ^ a, Fte — Fxi ^ a 
Ftm, flat A Fte, flat A 

As we will show in the next section, this linear decrease can be understood in terms of geometric optics. 



IV. COMPARISON WITH PERTURBATION THEORY AND GEOMETRIC OPTICS 



The aim of this section is to compare the numerical results of the previous section to those which were obtained 
from perturbation theory in Refs. EM E^ for an uniaxially and sinusoidally corrugated surface. We will show that 
discrepancies in the results from the two approaches can be qualitatively understood in terms of classical ray optics, 
a concept which was introduced in Ref. |2g for to the computation of Casimir interactions. In the perturbational 
path integral approach, the logarithm of the partition function is expanded in powers of the height profile hi as 
InZ = lnZ|o + lnZ|i + lnZ|2 + . . . . The zero order term lnZ|o = ^qALH^ 3 is the result for flat planes. The first 
order correction vanishes, lnZ|i = 0, since hi is on spatial average zero, and the second order contribution reads 

^ = 2^5/ h ^-j[ / f ^(|x||-xf||)[M*i)-M«i)] 9 - (40) 

where K (|x|| — xj| |) denotes a response kernel which has contributions from both TM and TE modes and was obtained in 
Ref. [2^. The second term is only finite for a smooth profile hi(xi) since the kernel has a singularity ~ |x|| — x'| |~ 3 . Thus 
for a rectangular corrugation with J [hi(xi) — hi(xi + x[)] 2 ~ \x[\ for \x[\ < A/4the perturbative result diverges due 
to the presence of sharp edges in the surface profile. In contrast, for a sinusoidal profile with hi(xi) = acos(27rai/A) 
one has J x [hi(xi) — hi(xi +x[)] 2 ~ and the divergence of the kernel is compensated. For this reason, we compare 
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our numerical results for the rectangular corrugation to the perturbative results for a sinusoidal profile 25, 26]. This 
will allow us to study the influence of edges on the Casimir interaction. Perturbation theory yields for the total 
Casimir force of the sinusoidal geometry of FiglBJa) the result 



F = F fla 



(41) 



with the function parameter free G(u) — ^r[5G(u) — uG'ju)] where G(u) = Gtm(u) + Gte(m) has contributions from 
TM and TE modes; for the explicit form of G(u) see Ref. |2(1 For comparison with our numerical results the limits of 
small and large H/X are of particular interest. From an expansion of G(u) one obtains the asymptotic expressions 

p ! — t 77 for A < H 




-Pflat 5fi) + Mf -20 (^)" forA»ff 



(42) 

In both limits the results are valid only if a -C A. In the limit of small A/a there is a divergence ~ a/A in the 
perturbative result which reflects the above mentioned divergence in Eq. H4U|I for rectangular corrugations with vertical 
segments. This singularity does not appear in our numerical results of the previous section; it is a characteristic feature 
of perturbation theory. In the following comparison we consider only the case A ^ a. Eq. I|42|) suggests for large plate 
separations H 3> A a decay of the excess force from the corrugation ~ a/H and for small H/X a decay ~ (a/H) 2 . 
The scaling behavior is in agreement with our observations for the rectangular corrugation as demonstrated by Fig|3] 
However, the latter Figure also shows that for smaller A/a < 10 the scaling regime with a decay ~ (a/H) 2 does not 
exist. 

Next, we will compare the perturbative results of Eq. I|42l) with our numerical results for the deviation of the actual 
Casimir force from the proximity force approximation (PA), (F — -FpA)/-Fflat> where Fp& is the force obtained from the 
PA. This approximation does no distinguish between the two types of modes and thus for the rectangular corrugation 
one has Fpa — 2Foo with Foo given by Eq. (|37|) . In general, for deformed surfaces the PA is ambiguous [2£| since the 
pairs of small parallel surface elements can be chosen to be parallel to either surface so that the local plate distance is 
measured either normal to Si or normal to S% as indicated by the arrows of Fig[5] We emphasize that this ambiguity 
does not arise for the rectangular corrugation. For smooth surfaces with finite curvature like a sinusoidal corrugation 
the PA result depends on the reference plate. If one measures the local distance perpendicular to the flat surface, as 
it is most common, one obtains for the Casimir energy per surface area 

£pa=-t [ dS£^[H - /n(zi)L (43) 
A Js 2 

but if the local distances are chosen perpendicular to the corrugated plate, one has 



1 f 

£pA,corr = ~7 I dS £fl a t 
A J Si 



(#-/ii(a:i))^l + (/ii(:ci)) s 



(44) 



where £fl a t(-ff) = — ^if -3 is the Casimir energy per surface area for two flat surfaces. For a sinusoidal corrugation 
the integrals over the surfaces can be computed perturbatively in a. This yields for large A the difference between the 
force F from perturbation theory [Ea. (|42Jl ] and the PA force, based on the flat and the corrugated plate, respectively, 



F — Fpa ^ 2 _ 20 W , *--fPA.co~ = r^!-20)(^) . (45) 



■Fflat \ 3 J \X' Fflat V 3 / \A 

The essential result is that the perturbatively obtained force approaches the PA approximation like (a/ A) 2 for large A 
which has to be compared to the a/ A decay seen in our numerical results for the rectangular corrugation, cf. Figs. El 
and Thus the deviation from the PA is stronger for the rectangular corrugation than for the sinusoidal profile, 
presumably due to sharp edges. Before we give a simple physical argument for the variation of the decay exponent let 
us compare the amplitudes in Eq. (|45[) . If we chose the PA to be based on the flat plate, the amplitude is negative, and 
the force Fpa is not a lower bound to the force at a fixed H/a, in disagreement with our observation for a rectangular 
corrugation. The corrugated surface based PA in contrast yields a positive amplitude. We expect that also for a 
sinusoidal corrugation the actual force is monotonous in A/a at fixed H/a, assuming its minimal value for A/a — > oo. 
The change of sign is just an other manifestation of the ambiguity in the proximity approximation. The observation 
that the actual Casimir force is located between the flat and the curved surface based PA was also made for a plane 
plate-sphere geometry recently |37^ . 
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FIG. 8: (color online) Typical paths of the proximity force approximation and the geometric optics approach for both sinusoidal 
(a) and rectangular corrugation (b) with A>a. Paths with arrows denote distances which are measured normal to one of the 
surfaces as used for the proximity force approximation. Paths without arrows denote the shortest surface connecting paths of 
length £(x) through a point x located in the gap between the plates. 



In order to understand the dependence of the exponent for the scaling towards the PA limit on the shape of the 
corrugation it is instructive to consider classical ray optics. Such an approach was recently applied to the calculation 
of Casimir interactions [2S|. Since this approach does not take into account diffraction it is limited to deformations 
where the radii of curvature are large compared to the smallest distance between the surfaces. But still, geometric 
optics allow for a better description of Casimir forces than the conventional proximity force approximation. By 
considering instead of all actual optical paths only the shortest paths, Jaffe and Scardicchio proposed an "optimal" 
proximity approximation for scalar field fluctuations subject to Dirichlet boundary conditions |2g. It can be also 
applied to electromagnetic fields. Consider a position x in the vacuum space between the plates, and denote by ^(x) 
the length of the shortest optical ray between the plates through that point. FigEI shows typical paths for the two 
types of corrugations we consider here. The Casimir energy in this optical approximation can then be written as 

P± = /rf 2 x,| T dx 3 f 3 (46) 

fflat J Ai(xr) A£ 4 (X||,IE3) 

where the integral runs over the total space between the surfaces. 

First, we apply this approach to the sinusoidal profile, see Fig|Sfa). For simplicity, we replace the sinusoidal profile 
by a piecewise linear profile, cf. Fig|5fa), which is a good approximation in the limit a <C A considered here. Then 
we have to determine ^(x) for each position between the plates for this simpler profile. Since the exact value for 
£(x) is difficult to evaluate, we consider the two cases where the position is close to one of the two surfaces and then 
assume a linear interpolation between the two length for £(x) at arbitrary x in the gap between the plates. If x is 
very close to the deformed surface Si the shortest paths is perpendicular to the flat surface 5*2. Contrary, if x is 
located close to the flat surface 5*2 the shortest ray is perpendicular to the deformed surface Si . With the so obtained 
approximative lengths ^(x) we obtain from Eq. I|4t)l) by expansion in a/H for the correction to the flat surface based 
proximity approximation the scaling behavior 



-Fflat 



Thus the optical approach nicely reproduces the correct scaling of the corrections to the proximity approximation at 
large A, in agreement with the perturbative result of Eq. I|45ll . 

In order to examine the role of edges for deviations from the proximity approximation, we apply the optical approach 
also to the rectangular corrugation of FiglHfb). For this geometry the shortest paths are easily identified. Except for 
positions located in an almost triangular shaped region (composed of the two shaded regions of FiglSfb) the paths 
are just perpendicular to both surfaces. Thus the deviation from the proximity approximation is caused by paths 
through points which are located inside the shaded region. These paths run either to corner C of the surface (larger 
region) or to the vertical surface segment (smaller region). For sufficiently large A the regions from adjacent edges do 
not overlap and can be treated independently. Furthermore, since the ratio of the area of the larger shaded region 
formed by the triangle ABC and the area of the smaller shaded region bounded by the vertical surface segment scales 
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like ~ (H/a) 2 one has to consider only the triangle ABC for the evaluation of Eq. (|46(l in the limit a/H <C 1. This 
gives 



opt 



Fflat V H A 



a a 



(48) 



This result is in agreement with the scaling behavior we have observed in our non-perturbative approach for the 
rectangular profile, see Figs El and and Eq. We conclude that the analysis of the shortest optical paths 

explains the observed dependence of the Casimir force on the surface shape close the proximity force limit A 3> H. 

Finally, we consider the ratio Ftm/Fte of the force contributions from TM and TE modes. In perturbation theory 
one obtains from the separate contributions of the two types of modes to the result of Eq. I|41|) the low a expansion 

which is valid if both H 3> A and X^$> a. Thus for sinusoidal corrugations the force has always larger contributions from 
TM modes at asymptotically large H, in contrast to our numerical results for rectangular corrugations, cf. Fig|5ja). 
We argued in the previous section that edges might cause the amplification of TE mode contributions. However, the 
convergence of the ratio to one for large H turns out be insensitive to the shape of the corrugations. Our numerical 
results agree perfectly over the full range of studied A/a with perturbation theory in that the ratio decays like a/H 
to one, see Fig|5fb). For small A/a — > the amplitude in no longer given by Eq. (|49() but saturates at a finite value 
which decreases with A since for A — > the reduced distance argument of section IIII Al implies equal contributions 
from both types of modes. 



V. CONCLUSIONS AND FURTHER APPLICATIONS 



In this paper we have developed a non-perturbative method to compute Casimir interactions in periodic geometries. 
This approach is based on a path integral quantization of the electromagnetic field subject to ideal metal boundary 
conditions. The so obtained effective action for the Casimir interaction is transformed to a representation which is 
adapted to periodic geometries and allows for an efficient numerical computation of the force between macroscopic 
objects. In particular, the approach allows us to compute the Casimir force between surfaces with strong periodic 
deformations and edges. For uniaxial deformations the electromagnetic field can be decomposed into two scalar fields 
which are subject to Dirichlet and Neumann boundary conditions, respectively. This enables us to study qualitative 
differences in the geometry dependence of the Casimir interaction for scalar fields with different boundary conditions. 
Applications of the latter case range from thermal fluctuations in superfluids to liquid crystals [2^ . l30| which can be 
described by a scalar field. Path integral quantization in the presence of boundaries has been previously applied to 
perturbative calculations of Casimir interactions between static and dynamic deformed manifolds in the context of 
both thermal 

mm 

and quantum fluctuations m,EL|2^,m] °f the confined field. However, all these computations 
were restricted to slightly deformed surfaces and edges were excluded. While a number of qualitative predictions 
of perturbation theory are confirmed by our approach even for strong deformations, we find novel non-perturbative 
effects which were unaccessible previously. 

As an explicit example, we calculated by the Casimir interaction between a flat and a rectangular corrugated plate 
with edges, including the case of large deformation amplitudes. Arbitrary periodic profiles can be treated by our 
approach as well by Fourier transforming the kernel of the effective action numerically and then applying the same 
technique we used here for the rectangular corrugation. We could confirm the perturbatively predicted existence 
of two different scaling regimes for the deformation induced part of the interaction as a function of the mean plate 
separation H . However, we also find that for small corrugation lengths only the large H scaling regime exists. We 
demonstrate by explicit calculations that in the limit of very small corrugation lengths the force can be obtained 
as the interaction of two flat surfaces with a reduced distance. At very large corrugation length and small H we 
find that the force approaches the result of the proximity force approximation. Our approach also allowed for a 
precise computation of the scaling of the force close to the limits of small and large corrugation length which provide 
an upper and lower bound, respectively, to the force. In both cases we find power law scaling with A/a, rendering 
corrections to proximity approximation in general large. The exponents of these power laws depend on the type 
of modes (transversal electric or magnetic) for small corrugation length. At large corrugation length we find an 
interesting dependence of the exponents on generic features of the corrugations. By comparison with perturbation 
theory for a sinusoidal corrugation we find that edges induce a slower decay towards the prediction of the proximity 
approximation as compared to smooth profiles. We could explain this generic behavior in terms of classical optical 
paths. 
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Our non-perturbative method can be applied to a number of other interesting situations. Since the path integral 
technique can be used in arbitrary dimensions of the embedding space and the surfaces, our method can be also 
used in this general cases. In this paper we focused on uniaxial deformations. Two directional corrugations can also 
be treated by our method by applying it to the full electromagnetic gauge field without splitting into TM and TE 
modes. The latter case could help to understand the possibility of repulsive forces since plates with two directional 
corrugations form at short distances cavities, i.e., geometrical shapes similar to a sphere for which a repulsive "force" 
is expected |l6j . At short plate separations, material properties become in general important for the interaction. 
These effects can be also described by a path integral approach with non-local boundary conditions (3^ , enabling the 
application of the methods developed here. For two corrugated surfaces, the existence of a lateral Casimir force has 
been predicted and computed by perturbative techniques [23, [2(| . It would be interesting to study the effect of strong 
corrugations and edges on the lateral Casimir effect by our method. For the dynamic Casimir effect the surfaces are 
dynamically deformed which leads for oscillations in time again to corrugated surfaces in Euclidean space but now 
along imaginary time. Our results thus imply different behavior at small and large frequencies. 



APPENDIX A: FOURIER TRANSFORM OF THE RECTANGULAR CORRUGATION MODEL 

We calculate the Fourier transformed matrices M. for Dirichlet and Neumann boundary conditions for a slightly 
more general geometry with two corrugated plates. Both plates are assumed to have a rectangular corrugation profile 
with the same wavelength A, but with different amplitudes ai and a^- This geometry is depicted in F'ig'EI an d the 
geometry of the system discussed in section III is obtained by simply setting the amplitude of the second plate to 
zero. The reason to perform this calculation here is that it is more transparent than the calculation which assumes 
one corrugated and one flat plate. In addition, we allow the plates to have a lateral displacement b. 

We start with the matrix for Dirichlet boundary conditions, cf. Eq. (|14af> . Performing first the Fourier transforma- 
tion with respect to xj_ = (xo, x 2 ), we have 

Jx.±Jy_ L Jx 1 Jy 1 

, (Al) 

= (27r) 2 (5 (2) (p ± +q ± ) / / / e *(Pi-Pi)*i+»(9i+Pi)»i 
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To evaluate this last expression analytically, it is necessary to find a simplified expression for the dependence of the 
second exponential term on x\ and y\. At this point, the use of piecewise constant profiles for the material plates 
becomes crucial: Since h a = ±a Q , for a = (3 we can write 

e -P~IM*i)-Mj/i)l = e - a ^[ch(a a p)+a- 2 h a (xi)h a (yi)sh(a a p)]. (A2) 
Similarly, for a ^ j3, we get 

e -p\h a ( Xl )- h ,(y l)+ H(S a2 -s (i2 )\ = e - P H [ ch ( aQ p)_(_i)« a -i/ la ( Xl ) sh ( aQ p)] [ c h(a^)-(-l)' 3 a^ 1 / l/3 (x 1 )sh(a^)]. (A3) 

To keep the notation short, we introduced p = \/p 2 L + pf- Now, we insert the Fourier series expression for h a given 
by 



7T — ' 9r* — 1 



n-1 



n— — 00 



into the rhs of Eqs. (JA2J1 and (|A3|) . Then, inserting those into Eq. IjAlJl . the remaining integrals over £1,2/1 and p[ 
can easily be performed. This yields the periodic formula 

00 

■M D (p,q) = (27r) 3 ,5^(p ± +q ± ) ^ S (pi + qi + 2nm/ A) Njy >m (q±, qi) (A5) 

m— — 00 

with the matrices 
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for m even, and 

, T , , / o c°i 2 (?±,gi) 

JVb,m(«±,«l) = D , . 

V C m,21 (?J-,?l) 

for m odd. The entries of the matrices are given as follows 



(-1) J 



(-1) 



E 

fe=— oc 
2i oo 

A;— — oo 



l 



-2a a q 2 k-l 



(m - 2k + l)(2k - 1) 



7 



(2fc-l)(5 (32 ^5 Q2 ) 



Q2k-1 

? — <l2k-\H 



(m-2k+l)(2k- 1) g 2 fc-i 



■ sh(a Q g 2 fe-i) sh(a/3Q2fc-i) ■ 



and 



Cm,a/3 (9J-)9l) 



mir 



(-1)"7 



i h(a ct9 )ch(a /39 ) + (-l)' 3 7 " 1 ^ 



• Bh(apq m ) ch(a 



(A7) 

(A8) 
(A9) 

(A10) 



where the phase factor 7 = e 2 " 6 was introduced. We note that the off diagonal entries a/3 and Q ^ implicitly 

depend on b through 7. Furthermore, g„ = \J q\ + {q% + 27rn/A) 2 was introduced, which implies q = qo- If a 2 = 0, 
the matrices A^.m have the symmetry Nu,m (<Z±j — Qi ) = Njy t — m (qj_, q%), and analogously for the Neumann matrices 
Nyi tm which we used in Sec. II. We remark that this symmetry is no longer valid for either type of boundary conditions 
if h 2 (xi) ^ h 2 (-xi). 

The matrix A^n for the Neumann boundary condition is obtained similarly as for the Dirichlet boundary condition. 
Evaluating first the Fourier transform of the orthogonal components as done in expression IjAlfl . the result is 

Mf (p,q) - (2n) 2 S^(p ±+q± ) 

e i PiXl+iqiyi { _ ir+ ( _ 9 2 3 + { h ' a ( Xx) + h ' ( yi )) dxidx3 - h' a ( Xl )h' p ( yi )d 2 Xl ) 



-ip[(xi-yi) £ 



-y/pl+P? \ X 3-V3\ 



X3=h a (x 1 )+H8 a2 
y3=hp(yi)+H5p 2 



iVpT+p? 

= (27r) 2 (5 (2) (p ± +q ± ) J I I e ^-v'xW+i{<ii+Pi 

Jx x J J/i Jp' 1 



(All) 



)yi 



(-1) 



Q+/3 



pI + p? 



VP±+Pi 



(d Xl -d Vl ) 



P? 



d Xl d yi 



We apply partial integration to obtain 



(_l)«+/3 



pi 



(Pi - 9i - Vi) 



Pi 2 



(Pl + P?) 



-{pi -p'i)(qi +p'i) 



(A12) 



i(pi-p' 1 )xi+i(gi+p' 1 )j/i p -\/p 2 ± +p'i \h c ,(x 1 )-hff(y 1 )+H{S c ,2-S l32 )\ 



Xi Jyi 



This expression will be treated analogously to the case of the matrix for the Dirichlet boundary condition, cf. 
Eq. IjAljl . It differs from the Dirichlet kernel by the additional p[ dependent term. This yields again Eq. I|A5|) . but 
now with -/Vd,™ substituted by the Neumann matrices N^ jTri , which are given by 



l m Bl , 21 (q±,qi) l m A*(q ± , qi ) 



. , . , § e qH ch(aig) ch(a 2 q) 



■2(l + e ^ 2ai<? ) 



JmO 



2 e -« if ch(aig)ch(a 2 ?) -|(l + e - 2a2 «) 



(A13) 
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H 



a 2 



r s 2 



*1 



Si 



FIG. 9: Two rectangular corrugated plates with the same wavelength A but different amplitudes ai and 0,2 and a lateral shift 
of b. The plates are translationally invariant along the X2 direction. 



for m even, and 



iV N , m {q±, qi) 



C£ tl2 (q ± , qi ) 
C^,2i(9±.9i) 



for m odd. The entries are now given by 

(-1)- 



A m , a {q±,qi) 



■ E 

k— — oo 



(m-2fc+l)(2fc-l) g| fc _! 



1 — e _2a c'J2fc-l 

~3 



(?_L,9l) , 



(A14) 



(A15) 



= 2 



(-1) 



7T 2 ,2^ 



7 



(2fe-l)(5g 2 -5a2) e -Q2k-lH 



k— — oo 



(ro-2fc + l)(2fc-l) g| fe _ i 



sh(a Q <7 2 fe-i) sh(a /3 g ; 2/ s -i) <j>mk{q±,qi) , (A16) 



and 



Cm, a/3 (<7-L,9l) 



m7r 



(-l) Q 7 m5a2 e- 9H + ^P^) sh(a a9 ) ch( a/3 g) 



+ (-l)^ 7 ^e-^ ff (g m - 



27TTO (/I + 27T77j/A\ 

A q m J 



(A17) 



sh(apq m ) ch(a a q m ) 



using the function 

<j>mk{q±,qi) = qi (qi + ^) (qi + y (2fc - 1)) + 2<?i ( 9l + ™) ^ + y (2fc - 1)) + q\ . (A18) 

As in the case of the Dirichlet matrices, the off diagonal elements depend on b via the phase factor 7 = e 2 ^ ±b . The 
matrices of the previous discussion of the rectangular corrugation model are now simply recovered by performing the 
limit a 2 — > and by denning a = a\. 



APPENDIX B: THE LIMIT OF SMALL A FOR THE MATRICES N m 

In this section, the limit A — > of the matrices N m (q±, qi) for the rectangular corrugation model of section III will 
be performed (cf. Appendix A for a — a\, a 2 — and A — ► 0). These matrices depend on the shift of the argument 
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qi relative to 2im/\ which requires a separate treatment of various cases. Considering this, for the Dirichlet case we 
find the simplified expressions 



A D , (q±,qi + 27m/A) A =° < 



e- Aaq + l e~ 



4q 



2q 



■ ch(aq) 



2q 



■ ch(aq) 



1 e~ 



2q 



e 

A 

47r|n| 



e 

4tt\u\ 



for n = 



for n odd • 



for n even 



(Bl) 



We have introduced a small quantity e, which is needed in order to have a non singular matrix Bki- However, at 
the end we can safely take e — > in the final expression for the Casimir force. As A — > 0, this quantity vanishes as 
e ~ Aexp(— 2irn(H — a)/X). The other matrices for m ^ are given by 



A D ,m (q±,Qi + 27m/ A) A =° < 



(-1)"^ e Z ^ sh( ) 

o o 



7rm 

o o\ 

/ 



■ sh(ag) 



for n = 



for n = —m 



for n £ {— m, 0} 



(B2) 



for m odd, and 



Nu.m (q± , qi + 2tto/ A) A =° < 



(-l)T e" 2 "" 
7T 2 n(m+n) q 





± 








for n odd 



for n even 



(B3) 



for even m ^ 0. Analogously, for the Neumann matrices, we find 



' / -f(e- 2a 9 + l) |e-^ ch ( ag ) \ 
{ ^e-^ch(aq) -§ J 



for n = 



A N ,o (q±,qi + 27m/ A) A =° < 



4(-l)"- 1 g 2 / -2ag 



^7 



7r|n| 



7r|nJ 
A 



for n odd 



for n even 



(B4) 
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and 



2(~l)^~ qi p -qH 



sh(aq) 



for n = 



Nist,m (q± , qi + 27rn/ A) 



A-+0 



2(-l) ^ gi -off 



e-^sh(ag) 








for n = — m 



for n $ {—to, 0} 



(B5) 



for m odd, and 



A^n,™ qi + 2wn/ A) A =° < 



A 2 9 3 







2n(n+Ti) 



— C m (n) 




±^C m (n) 



for nodd 



for neven, n ^ {— m, 0} 



for n € {— m, 0} 



(B6) 



for even m ^ 0. Here, the asymptotic behavior of e for A — > is e ~ A 1 exp(— 2im(H ± a)/A). The constant is given 
by C m (n) = (— l) m / 2 X];=-oo[(^^ — -OC^ — 1 — m) 1 2Z — 1 + n\ and the prime at the summation sign indicates that 
I ^ (l-n)/2 if n is odd. 



ACKNOWLEDGMENTS 



This work was supported by the Deutsche Forschungsgemeinschaft through Emmy Noether grant No. EM70/2-2. 



H. B. G. Casimir, Proc. K. Ned. Akad. Wet. 51, 793 (1948). 

P. W. Milloni, The Quantum Vacuum (Acad. Press, 1993). 

M. Bordag, U. Mohideen, and V. M. Mostepanenko, Phys. Rep. 353, 1 (2001). 

K. A. Milton, The Casimir effect: Physical Manifestations of Zero-Point Energy (World Scientific, 2001). 
M. Kardar and R. Golestanian, Rev. Mod. Phys. 71, 1233 (1999). 

J. Israelachvili, Intermolecular and Surface Forces (Academic Press, San Diego, 1992). 
R. Garcia and M. H. W. Chan, Phys. Rev. Lett. 88, 086101 (2002). 

A. A. Bytsenko, G. Cognola, L. Vanzo, and S. Zerbini, Phys. Rep. 266, 1 (1996). 
K. A. Milton, Phys. Rev. D 22, 1441 (1980). 

K. A. Milton, Phys. Rev. D 22, 1444 (1980). 

S. K. Lamoreaux, Phys. Rev. Lett. 78, 5 (1997). 

U. Mohideen and A. Roy, Phys. Rev. Lett. 81, 4549 (1998). 

H. B. Chan, V. A. Aksyuk, R. N. Kleiman, D. J. Bishop, and F. Capasso, Science 291, 1941 (2001). 
G. Bressi, G. Carugno, R. Onofrio, and G. Ruoso, Phys. Rev. Lett. 88, 041804 (2002). 

B. Derjaguin, Kolloid Z. 69, 155 (1934). 

F. Boyer, Phys. Rev. A 9, 2078 (1974). 

D. Iannuzzi, I. Gelfand, M. Lisanti, and F. Capasso, preprint quant-ph/03 12043 

G. J. Maclay, Phys. Rev. A 61, 052110 (2000). 

F. M. Serry, D. Walliser, and G. J. Maclay, J. Microelectromech. Syst. 4, 193 (1995). 
A. Roy and U. Mohideen, Phys. Rev. Lett. 82, 4380 (1999). 

F. Chen, U. Mohideen, G. L. Klimchitskaya, and V. M. Mostepanenko, Phys. Rev. Lett. 88, 101801 (2002). 

G. L. Klimchitskaya, S. I. Zanette, and O. A. Caride, Phys. Rev. A 63, 014101 (2001). 
R. Golestanian and M. Kardar, Phys. Rev. Lett. 78, 3421 (1998). 



21 



[24] R. Golestanian and M. Kardar, Phys. Rev. A 58, 1713 (1998). 

[25] T. Emig, A. Hanke, R. Golestanian, and M. Kardar, Phys. Rev. Lett. 87, 260402 (2001). 

[26] T. Emig, A. Hanke, R. Golestanian, and M. Kardar, Phys. Rev. A 67, 022114 (2003). 

[27] R. Balian and B. Duplantier, Ann. Phys. (N.Y.) 112, 165 (1978). 

[28] R. L. Jaffe and A. Scardicchio, preprint quant-ph/0310194 

[29] H. Li and M. Kardar, Phys. Rev. Lett. 67, 3257 (1991). 

[30] H. Li and M. Kardar, Phys. Rev. A 46, 6490 (1992). 

[31] T. Emig, Europhys. Lett. 62, 466 (2003). 

[32] M. E. Peskin and D. V. Schroeder, An Introduction to Quantum Field Theory (Addison- Wesley, Reading, MA, 1995). 
[33] J. D. Jackson, Classical Electrodynamics (John Wiley & Sons, Inc., 1967). 

[34] S. K. Karepanov, M. Y. Novikov, and A. S. Sorin, Nuovo Cimento Soc. Ital. Fis., B 100, 411 (1987). 
[35] G. L. Klimchitskaya and Y. V. Pavlov, Int. J. Mod. Phys. A 11, 3723 (1996). 
[36] V. M. Mostepanenko and N. N. Trunov, Sov. Phys. Usp. 31, 965 (1988). 
[37] H. Gies, K. Langfeld, and L. Moyaerts, JHEP 06, 018 (2003). 
[38] R. Biischer and T. Emig, preprint cond-mat/0308412 



